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ABSTRACT 


Work performed during the first twelve months of the project duration 
for the NASA Grant (NAG-1-861) is reported here. An analytical and compu- 
tational study of opposed jet diffusion flame for the purpose of under- 
standing the effects of contaminants in the reactants and thermal diffusion 
of light species on extinction and reignition of diffusion flames is in 
progress. The methodologies that have been attempted so far are described. 
Results using a simple, one-step reaction for the hydrogen-air counterflow 
diffusion flame are presented. These results show the correct trends in 
the profiles of chemical species and temperature. The extinction limit can 
be clearly seen in the plot of temperature vs. Damkohler number. 


INTRODUCTION 

The laminar, counterflow diffusion flame has become a standard tool for 
the investigation of stability and extinction characteristics of diffusion 
flames. In one configuration, the counterflow diffusion flame is obtained 
in the laboratory by means of flow from two coaxial jets -- fuel from one 
side and air from the other -- impinging upon each other and forming a 
stagnation point flow. In a different configuration, air flows towards a 
porous cylinder from whose surface fuel is injected. These two streams 
form a stagnation point flow at some distance from the cylinder surface. 

The two configurations described above have their advantages and disadvan- 
tages. From the analysis point of view, both are described by the classi- 
cal boundary -value problem for stagnation point flow. They differ only in 
the boundary conditions. The present work describes a procedure for ana- 
lyzing problems of the boundary -value type. Even though the method is 
quite general, attention will be focused here on the counterflow difusion 
flame . 
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Extinction of the counterflow diffusion flame has been studied exten- 


sively because of its importance in predicting blowoff conditions in a 
practical jet engine. This aspect is becoming increasingly important as 
new engine concepts such as the supersonic combustion ramjet (SCRAMJET) 
engine are being developed. A non-dimensional parameter called a Damkohler 
number is convenient for the study of extinction characteristics of non- 
premixed flames. We difine a Damkohler number, , as the ratio of the 
flow time to chemical reaction time. Therefore, smaller Damkohler numbers 
represent slower reactions and larger Damkohler numbers indicate faster 
reactions in relation to the flow time. Extinction and reignition phe- 
nomena have been studied with respect to a plot of the flame temperature as 
a function of the Damkohler number. An S-shaped curve have been obtained 
for a limited number of cases (Figure 1) . 



Figure 1. Graphical Representation of Flame Limits 
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Stability of the flame can be studied with respect to this S- shaped 
curve. The two vertical tangents at points 1 and 2 represent the extinc- 
tion and reignition points. The branch A of the curve represents a stable 
flame; the branch B represents unstable bumig conditions and branch C 
represents another stable flame. Even though this S- shaped curve have been 
obtained analytically by asymptotic analysis, the existence of it has not 
been established by more detailed solution of the governing equations; 
nor has it been experimentally demonstrated in a universal manner. The 
motivation for the present study is to conduct a more thorough analysis 
than has been done before and answer some of the interesting questions that 
have been raised by the experimental work of Pellet et al[l] . Some of the 
interesting observations of Pellett and others are decribed in the follow- 
ing paragraph . 

An argon-bathed tubular opposed jet burner (OJB) was developed and 
applied to chemically- limited counterflow diffusion flame, formed by small, 
( 2.7 and 5 mm i.d. ) laminar jets hydrogen diluted with nitrogen and air 
at ambient conditions. Data from pure air as well as air mixed with 
nitrogen were characterized and used to evaluate similar data for nitric 
oxide contamination. Two phenomena called 'blowoff' and 'restore' were 
observed by varying test parameters such as fuel dilution ratio and jet 
velocity. 'Blowoff' was a sudden breaking of the disk- shaped flame to form 
a stable torus shape; 'restore' marked the restoration of the flame to the 
center by radial flame propagation. 'Restore' is a newly discovered flame 
limit. Both flame limits were sensitive to reactant composition. Blowoff 
is a measure of flame reactivity and it may be related to a Damkohler num- 
ber defined as the ratio of flow time to reaction time. 'Restore', on the 
other hand, appears to be more complex and cannot be explained by flame 
reactivity alone. The inner edge of the torus -shaped flame is approxi- 


3 



mately the same diameter as the tubes from which the reactants flow. It is 
believed that 'restore' is a flame propagation phenomenon: when the fluid 
velocity at the inner edge of the torus -shaped flame is less than the flame 
speed, the flame propagates to the central region marking the phenomenon of 
restore. However, no conclusive evidence is avaialble to back this reason- 
ing. The question that still remains to be answered is whether mixing by 
diffusion is complete by the time the gases travel outward a distance of 
one tube diameter. In any case, there is reason to believe that the coun- 
terflow flame is a diffusion flame only in the neighborhood of the stagna- 
tion point; away from the stagnation point it exhibits some characteristics 
of premixed flames. Therefore, a simple diffusion flame theory may not be 
adequate to explain all the phenomena observed in the laboratory. 


PRESENT WORK 


The present work consists of the solution of the equations governing 
the fluid flow, conservation of energy and finite rate chemistry for the 
hydrogen-air counterflow diffusion flame. The equations are reduced to a 
set of ordinary differential equations by making the classical boundary 
layer assumptions; they are, therefore, valid only for the stagnation point 
flow. These equations are of the boundary-value type. The combustion 
model in the preliminary calculations assumes a one -step forward overall 
gas -phase chemical reaction of second order occur ing in the axisymmetric or 
two-dimensional stagnation point flow at a solid surface or formed by two 
opposed jets. The resulting set of equations is as shown below. 


(c f„), + + oh) f 

{ k ) f ^ ~ 0 
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The chemical reaction rate for the one -step second order reaction is 
given by 

** = z, y 0 C- * /*t) 

in which w is mass of fuel consumed per unit volume per unit time. 

Define the first Damkohler number as 

T> = ^ t /*- 

Introducing Eq. (5) into Eq. (3), the species conservation equation becomes 

(- L y Fn ) + t Vf, + J>YrYo (-Z =0 fsD 

/V 1 


Boundary Conditions 

The appropriate boundary conditions have to be applied depending on the 
geometry of the problem. For the opposed jet problem the boundary con- 


dtions are given as 
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f = o 
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For the stagnation point flow over a surface one can write the boundary 
conditions as 


f IcO = o 
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NUMERICAL SOLUTION 


Several different numerical approaches have been attempted; some of 
them proved to be more suitable for the present problem than the others. 

The main difficulty in the solution procedure stems from the fact that the 
equations become stiff as the number of chemical reactions increase. Stiff- 
ness is defined as the condition where the eigenvalues of the Jacobian Mat- 
rix of the locally linearized representation of the equation set have nega- 
tive real parts that fall over a large range; eigenvalues differing by 5 or 
6 orders of magnitude are common in chemically reacting flow problems. 
Stiffness limits the choice of step size that can be used for the numerical 
integration of the problem; as a result computation time becomes prohibi- 
tive. Several techniques are in use to circumvent the difficulties associ- 
ated with the solution of stiff equations. Multistep methods such as used 
in the widely used routines like LSODE [2] are useful for solving initial 
value problems in ordinary differential equations. No such well proven 
methods are available for boundary value problems in ordinary differntial 
equations or for solving stiff partial differential equations. The present 
work is based, therefore, on methods chosen by intuition rather than any 
well established criteria. 

Numerical algorithms for the governing equations described earlier have 
been developed and they have been implemented on the NASA Langley computer 
system as well as the computer system at the University of Missouri-Rolla. 
The following are brief descriptions of the algorithms. 

1. Collocation Method: This is a routine called COLSYS [3] available from 

the Langley Math Library. The routine can be used to solve boundary- value 
problems governed by a system of ordinary differential equations. 

2. Fully Implicit Finite Difference Method: Recognizing the stiff nature 

of the differential equations because of the presence of the chemical 
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source term, a fully implicit algorithm was written to avoid the limitation 
on step size for stable solution. 

3. Patankar and Spalding Method: This method uses the conservative form 

of the equations. The source term is treated according to the rules for- 
mulated by Patankar [4] . 

Work in Progress 

The method that has been used for the results given in this report is a 
combination of COLSYS for the solution of the fluid dynamic equation and 
the fully implicit method for the solution of the energy equation and the 
species continuity equation. The COLSYS routine has been found to be well 
suited for the solution of the fluid dynamic equation which is third order. 
Also it is well recognized that, for this case, the fluid dynamic equation 
is only weakly coupled with the energy and the species continuity equa- 
tions. Therefore, the fluid dynamic equation is solved first with assumed 
initial profiles. With the values of the similarity variable f and its 
derivatives obtained from the solution by COLSYS, the energy equation and 
the species continuity equations are solved simultaneously by using the 
fully implicit solution procedure. The grid spacing in the eta direction 
and the time step are both variable. The grid spacing in the eta direction 
is chosen such that they are closely spaced in the region where the temper- 
ature and the species mass fractions have high gradients. The time step is 
varied in an ad hoc manner so as to reduce computation time and to preserve 
time accuracy of the solution. 
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PRELIMINARY RESULTS 


Some preliminary results of the calculations are shown in Figures 1 -5. 
These are preliminary results and therefore, their numerical values are not 
to be stressed. It was found that the results depend on factors such as 
the grid spacing. A grid independent solution should, therefore, be 
obtained by further refining the grid. A course grid is presently used to 
reduce computer run time. Variable time step was used to decrease the time 
required to reach asymptotic values. The results shown are the values 
obtained at the end of 3000 time steps. It was observed that slight oscil- 
lations of the solutions persist even at this stage; however, these were 
observed to be too small to have any significant effect on the solution. 

The CPU time required for these runs is approximately 50 minutes on the IBM 
4381 mainframe computer. 

Figures 1, 2 and 3 are for the hydrogen- air case. Figure 1 shows the 
temperature profile for one value of the velocity gradient. The peak tem- 
perature of 3300 K is higher than the adiabatic flame temperature because 
these values are based on constant specific heat assumption. Subsequent 
runs with variable specific heats have shown that the peak temperature is 
lower than the adiabatic flame temperature. Two other temperature profiles 
corresponding to extinction conditions are also shown in the plot. 

Figures 2 shows the corresponding mass fraction profiles. These pro- 
files agree in shape with published data. It can be seen that hydrogen and 
oxygen concentrations go to zero at the flame zone. Water exists in sig- 
nificant amounts over a wider region to either side of the flame zone. 

Figure 3 shows the extinction conditions. As the Damkohler number 
decreases the maximum temperature decreases; the vertical portion of this 
curve may be thought of as the extinction condition. Again, it must be 
emphasized that the numerical values are quite different from previously 
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published numbers because the values of the frequency factor and the acti- 
vation energy used in the calculations were not realistic; subsequent cal- 
culations show extinction strain rates much closer to experimentally 
observed values. 

Figures 4 and 5 show results for a carbon monoxide -air flame. The 
trends are similar to the hydrogen-air flame. As expected, the peak tem- 
perature is lower for this case. 


WORK IN PROGRESS 

The following is a list of cases that are being run presently. 

1. Opposed jet case. 

2. Effect of activation energy and frequency factor on extinction. 

3. Effect of time step size on the solution. 

4. Effect of grid spacing in the eta direction. 

5. Runs to compare the strain rate vs. fuel mole fraction for extinc- 
tion curve of Pellet et al. 
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Species Mass Fraction Distribution 
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